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Abstract 

To resolve infrared problems with the de Sitter invariant vacuum, we 
argue that the history of the de Sitter phase is crucial. We illustrate how 
either (1) the diagonalization of the Hamiltonian for long-wavelength modes 
or (2) an explicit modification of the metric in the distant past leads to nat- 
ural infrared cutoffs. The former case resembles a bosonic superconductor 
in which graviton-pairing occurs between non-adiabatic modes. While the 
dynamical equations respect de Sitter symmetry, the vacuum is not de Sit- 
ter invariant because of the introduction of an initial condition at a finite 
time. The implications for the one-loop stress tensor and the production of 
particles are also discussed. 
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1 Introduction 



One of the most perplexing problems of cosmology concerns the magnitude of the cosmo- 
logical constant,! for which a stringent observational upper limit exists. Although super- 
symmetry requires the cosmological constant to vanish, supersymmetry is either explicitly 
or spontaneously broken in nature, i As a result, the low-energy effective field theory is 
not supersymmetric, and the natural scale of the cosmological term in the gravitational ac- 
tion should be no smaller than the magnitude of the vacuum energy associated with the 
quanta of matter and gravitons themselves. As the universe cools, such a cosmological term 
would eventually predominate over radiation and other forms of matter. Indeed, the most 
popular cosmology today, the inflationary scenario (for a review see the book [Q), assumes 
that our universe passed through an era in which the cosmological term dominated, and it 
is a total mystery why, when this phase terminated, we should be left in a universe with a 
tiny or vanishing vacuum energy. We do not advance a solution to the cosmological con- 
stant problem here, but address the simpler but nevertheless vexing question of what is the 
correct vacuum state during the inflationary era. The answer to this question is a prelude 
to addressing consistently the issue of whether such a phase is intrinsically unstable due to 
gravitational fluctuations, as some have suggested. 0, [5], |6], |7|, |J It also is crucial for under- 
standing the quantum field theory of a free, massless, minimally coupled scalar field in a 
de Sitter background. [§J 

The general form for the action may be written as the sum of gravitational and matter 
terms: 



3 Spontaneously broken global SUSY demands a nonzero vacuum energy in the broken phase. In super- 
gravity a nonvanishing p vac is not obligatory and one can adjust the parameters so that it vanishes when 
the symmetry is broken but the required fine-tuning is unnatural and should be accurate with the precision 





2 For a review, see e.g. Ref. g|. 



of about 10 
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The gravitational action S g with a cosmological constant A is 

S g = ^ I d 4 x^(R -2A), (2) 

with k, 2 = lQ-K/m 2 Pl . The corresponding field equations are 

R»v-\g^(R-2K) = yT^. (3) 

where is the stress tensor of matter. It is presumed that the cosmological term, which may 
actually arise in part from vacuum fluctuations of particles, dominates the stress tensor of 
these particles, so we may neglect the matter action and T^ u for the moment, de Sitter space 
is the maximally symmetric solution of the resulting equations. In conformal coordinates 
appropriate to the spatially flat sections, the metric takes the form 

g»u = a 2 (r)v^u, (4) 

where rj^ is the usual Minkowski space metric, a(r) = —1/Hr, with H = \jA/3. 

The physical fluctuations in the metric, that is to say, the gravitons, in the neglect of 
self-interactions, obey the same field equation as the free, massless, minimally coupled scalar 
field. Although our ultimate interest is in the purely gravitational case, for simplicity, we 
will consider the quantum field theory of a massless, minimally coupled scalar field in the 
de Sitter background. This field obeys the same equation of motion as certain modes of 
the graviton field, including the physical modes, but avoids the intricacies of tensor algebra, 
gauge fixing, unphysical modes and ghosts that complicate the treatment of the graviton 
field itself. In this case, it is well-known that a de Sitter invariant vacuum state does not 



exist. |T0|, U To establish notation, we will review the situation briefly. The action for the 
scalar field is simply 

S m = I tfxyfg^* . (5) 



In the conformal coordinates introduced above, the field equation takes the form 

2d T a „ , , „ 2 



j; + _j-a r ^ = ^--a r = o. (6) 



The solution of this equation with definite momentum k is 



(f) k (r, x) = fc (r) exp(ik ■ x); 



with 0fc(r) = H ( r — — ) exp(— ikr) 



(7) 



The corresponding field operator may be written as 

d 3 k 



(27r)V2fc 



a^ fe (r)e^+a/0(r)*e-^ 



(8) 



t 



with a fc and a fc satisfying the commutation relations 



(27r) 3 5 (3) (£-g) 



(9) 



All other commutators are zero. The natural, de Sitter invariant no-particle state |0, in) 
that suggests itself is that state annihilated by the af 



|0, in) = 0. 



(10) 



Unfortunately, although this would appear to be a perfectly acceptable, adiabatic in-state, 
this choice of vacuum leads to infrared divergences, illustrating the difficulties first pointed 
out in Ref. [ID]. For example, the Feynman propagator in this state is ill-defined because 
of an infrared divergence. The source of the problem is the i/k term in the wave function, 
which is a feature not present in Minkowski space (or for the conformally coupled massless 
scalar in curved space.) The propagator takes the formf7| 



Gf{x, x') 



4tt 2 



-, ttt hi ur — x 

(x-x') 2 2 v ; 



(11) 



The scale of the argument of the log is not specified, reflecting the absence of an infrared 
cutoff. This propagator grows for large spacelike separations; simply inserting a scale into the 
logarithm would not change that. In an acceptable vacuum state, such correlations should 
fall to zero at large distances (so-called cluster decomposition.) It is noteworthy that the 
commutator of the fields does not display such infrared problems, a feature that we have 



exploited elsewhere. fllT|1 This emphasizes that this is a uniquely quantum effect, without a 



classical counterpart. One may object that the Feynman propagator is irrelevant, because of 
complications defining the out-state in these coordinates. But perturbation theory requires 
the computation of other interesting quantities such as the time-ordered product between 
in-vacuua or, similarly, the anticommutator G\ defined by 

Gi(t, t', x — x') = (0, m|{0(r, x), 0(r', x')}|0, in). (12) 

These, like the Feynman propagator, suffer infrared divergences. This infrared problem is 
associated with the behavior at large distances and is not peculiar to large times.B 

A useful discussion of this problem has been given by Allen and Folacci,[12| who concen- 



trate on a different foliation of the de Sitter manifold, viz., spatially closed coordinates. The 
advantage of these as opposed to flat (or open) coordinates is that the modes are discrete. 
Thus, the zero mode may be separated and modified, and they introduced modified vac- 



uua that were suggested by a construction of DeWitt [|13|| in which the zero mode is treated 
differently. However, such a construction does not resolve infrared problems for the flat or 
open coordinates, in which the momentum is continuous. Modification of the treatment of 
the zero mode alone will not change the cluster properties of the correlation functions. Con- 
sequently, people resort to the introduction of an ad hoc infrared cutoff (such as \k\ > k ) 
whose physical origin is uncertain. Although the stress tensor is not infrared divergent, there 
is little hope that all observables will be cutoff- independent. Even for the closed coordinates, 



it is not so clear that the so-called Hadamard vacuua discussed in Ref. JL2] are physically 
selected. 

So the question remains how this infrared problem is resolved physically. The existence of 
a physically reasonable propagator (i.e., two-point correlation functions) is a prerequisite to 
formulating a sensible perturbation expansion. Only then can one begin to ask, for example, 
whether de Sitter space is stable to quantum fluctuations. 



4 This divergence is, however, the origin of the instability discussed in Ref. 0. Other infrared problems 
at large times are the basis for the instabilities discussed in Refs. || @, §|. 
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The outline of the remainder of this paper is as follows: in the next section, we discuss 
the treatment of the nonadiabatic modes. In Section 3, we evaluate the change in the stress 
tensor for our new vacuum state. In Section 4, we consider a different method of treating 
the history, by modifying the metric in the distant past, and using a WKB expansion for an 
evaluation of the effects. In Section 5, we consider yet another modification of the metric, one 
which lends itself to analytic solution. In Section 6, we summarize our results and comment 
on their significance. 

2 Nonadiabatic modes 

Even though the stress tensor is known to be free of infrared problems, an examination of the 
Hamiltonian does in fact provide some insight into the nature of the problem. Because the 
metric is explicitly time-dependent, momentum but not energy is conserved. Nevertheless, 
even though the Hamiltonian is not time-independent, it still governs the time-evolution 
of the field. It is useful to construct the Hamiltonian corresponding to the matter action 
Eq. d). We find 



H{T) = / ^ k 



Sk ■ (,. ,. tx, t„ \ '1. „ >l „t„t 



2 \ k k k k J 2 k ~ k 2 k ~ k 



(13) 



where 



The coefficients r^, Sk depend only on the magnitude k = \k\. Notice that they are in fact 
functions of the product kr only. From this expression, we see that in general, the Hamilto- 
nian causes transitions increasing or decreasing the number of quanta by two. In particular, 
the Hamiltonian makes transitions between the de Sitter invariant vacuum |0,m) defined in 
Eq. flTUp and the two-particle states 

\k,-k)=a3a^\0,in). (15) 
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Note that the this form of the Hamiltonian is precisely that of a BCS superconductor in 
quadratic approximation,! except that the quanta here are bosons rather than electrons. 
Nevertheless, one is tempted to diagonalize the Hamiltonian at some particular time by a 
Bogoliubov transformation |14, [0| and to try to define no-particle states with respect to the 



transformed basis. This naive expectation is not correct, and the vacuum state would be 
extremely sensitive to the time chosen. Such a procedure has been a highly disreputable 
thing to suggest for many years, [[U], |17| especially because it leads to an infinite density 



of particles at any later time. In any case, the off-diagonal Sk vanishes in the distant past 
(r — > — oo), so one would be inclined to choose the in-state at r = — oo to be precisely |0, in), 
especially since the rate of change of the metric in conformal coordinates is so slow in the 
distant past, and the Hamiltonian approaches a finite, time-independent limit. 

This line of reasoning is perfectly sensible and would be the end of the story were it not 
for the fact that the theory remains infrared divergent in that state. Moreover, in our present 
thinking about cosmology, in which the inflationary phase arises either after cooling from a 
hot initial phase or from some sort of phase transition at very early times, in fact it is an 
idealization to presume one is in the de Sitter phase at all times. Let us consider how our 
description might change as a result of the supposition that the idealized Lagrangian that 
we have been discussing is only valid for times exceeding some early time r . Our concern is 
with what the state of the universe must be. At first, one would think that the no-particle 
state will correspond to the usual adiabatic vacuum.! In this construction, at any given time, 
the modes of the field may be divided into high- and low-momentum modes, the transition 
being set by the time rate of change of the scale factor. In conformal time, this is 

log(a(r)) = -. (16) 

GST T 

Since we may not entertain times earlier than tq, we may think of those momenta with 
A; I tq I < 1 as non-adiabatic and those with k\ro\ > 1 as adiabatic. For the adiabatic modes, 



5 See any modern text on superconductivity. [fTi] 



3 The adiabatic vacuum was developed by Parker[18] and others and is reviewed in Ref. |jl 



we will assume that the standard construction prevails, so that, to a good approximation, 
the correct vacuum state will be like |0, in) in that it will be annihilated by the at, but only 
for k > 1 / 1 tq | . 

However, the nonadiabatic modes require further discussion. It is reasonable to ask 
whether an observer might design an experiment to probe these low frequency or very long 
wavelength modes. In fact, because of the horizon in de Sitter space, they are impossible to 
detect. We allude to the well-known fact is that, no matter how long one waits, an observer at 
the origin can only receive signals that originate later than r at comoving coordinates \x\ < 
| tq | - Correspondingly, by the uncertainty principle, a detector could not be built to probe 
(conformal) momenta \k\ < 1/|t |, since no signal can be communicated over coordinates 
greater than |r |. Consider in particular the so-called two-particle states defined in Eq. flHf) . 
Normally, we think of them as consisting of two distinct quanta moving in opposite directions. 
However, for frequencies k < 1/|t |, it would be impossible to make such a determination. 
In fact, such nonadiabatic modes are essentially indistinguishable from the no-particle state 
|0, in) itself. Thus, we have a situation in which the Hamiltonian causes transitions between 
physically indistinguishable states, the classic situation discussed by Lee & Nauenberg. Hl9| 
They argued quite generally that one must use techniques of degenerate perturbation theory 
if one is to avoid encountering infrared divergences. In particular, they showed that one must 
either sum over degenerate states in forming observables or diagonalize the Hamiltonian 
in the space of degenerate states. In the present context, this suggests that, while the 
definition of the vacuum |0,m) given in Eq. fllOl) may be correct for the adiabatic modes, 
one should modify the vacuum for the nonadiabatic modes by reference to a basis in which 
the Hamiltonian is diagonal. 

In an interacting field theory, diagonalizing the Hamiltonian is a formidable task, but 
since we are working here with a quadratic Hamiltonian, it can be carried out explicitly by 
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means of a Bogoliubov transformation 

h = ak H + & a Jk> ( 17 ) 

with coefficients subject to the constraints 

|a fe | 2 - \/3 k \ 2 = 1. (18) 

ctfc and depend only on the frequency k — \k\ and not the direction.B 

This discussion suggests the following improved approximation to the vacuum state of 
de Sitter space. For an arbitrary Bogoliubov transformation, Eq. (|T7|), one may define a 
vacuum state |r ,m) by 

b % \T ,in) = 0. (19) 

The state is defined by the prescription one gives for the coefficients a k and (3 k - As suggested 
earlier, we choose the usual prescription = 1, (3k — for adiabatic modes k > 1/|t |. Thus, 
for these modes, this vacuum is the same as the de Sitter vacuum. However, for the non- 
adiabatic modes, k < 1 / 1 tq | , we choose these coefficients such that they diagonalize the 
Hamiltonian. (We shall come back to discuss the transition between the two regimes later.) 
While the physical motivation based on the Lee-Nauenberg theorem seems quite compelling, 
it remains to be shown that this prescription is precisely of the sort described by Ford and 
Parker [|l^], i.e., it resolves the infrared problems associated, for example, with the two-point 



correlation functions. In particular, since the Hamiltonian varies with time, it may not be 
immediately obvious that diagonalizing the Hamiltonian at some initial time tq will remove 
the infrared divergences for all later times. 

As noted earlier, the Hamiltonian is explicitly a function of time, but it is the Hamiltonian 
H(tq) at the initial time tq that we wish to diagonalize. It is a simple exercise to obtain 
the form of the Bogoliubov transformation. We define an angle 8 k by the phase of Sk, i.e., 



7 We will display the coefficients ak , Pk below. 
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Sk = \sk \ exp(—i9k). Then we find 



— = Rkexp(—i9 k ), with R k = — — 

Oik ' \Sk\ 



\ 2 

A -1. (20) 

\ s k \ 



A useful form for the phase angle is 

6k = 2kr — 7 fc + 7r, where tan7 fe = 2kr . (21) 

Without loss of generality, we may choose the coefficient a k to be real, so that the phase 
— 9 k is associated with (3 k . Then, for the magnitudes, we find 

<*k= r^— 2 > IAI= r^—^ fc < 1/|to|. (22) 
1 — Ri, \/l — Rh 



The Hamiltonian then takes the form 



^)=/(03*K 6 U» ( 23 ) 



As remarked earlier the state so constructed resembles superconductivity, as the vacuum 
state involves pairing of quanta of equal and opposite momenta. It is important to recognize 
that the pairing of nonadiabatic modes does not correspond to condensation of the scalar 
field. The symmetry breaking arises from off-diagonal long-range order. Thus, in terms of the 
original quanta, there are extremely long-range correlations in the vacuum which resolve the 
infrared divergences. Note, however, that unlike ordinary superconductors, diagonalization 
does not result in a mass for the elementary excitations. Of course, the expectation value of 
the Hamiltonian H{tq) is lower in the state |r ,m) than in the original de Sitter invariant 
state |0,m). So we may think of this as spontaneous breaking of de Sitter symmetry, but it 
also exhibits aspects of explicit breaking, since it involves a parameter (r ) that is not part 
of the original Lagrangian. On the other hand, if we regard the gravitational action as being 
defined only on the interval of time To < r, then one may say it was present in the theory 
from the beginning. 
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In the transformed basis, the expansion of the field becomes 



d 3 k 



(27r) 3 V2A; 

where the modified wave functions are 



(24) 



(25) 



It is useful to have approximate expressions for these new wave functions for small momenta. 
In the nonadiabatic regime, k\r \ <C 1, we then find R k ~ 1 — 2(kr ) 2 , so that 

1 / 1 



Oik 

Ok 



2k\r 
7i + ^(kr ) 2 + . . . 



(26) 



In particular, a k + (3 k — > O(k\r \) as k — > 0, as required for removable of the infrared 
singularity. [K| To see this explicitly, we may use these results to obtain the behavior of the 
wavefunction, Eq. (p5"D, in the nonadiabatic regime. Since the conformal time is in the range 
T < t < 0, we may expand for all time to obtain 



u k (r) w iHr 



1 + - {{krf + (kT f) 



(27) 



Thus, the \ jk infrared behavior as k — > characteristic of the de Sitter wave function, Eq. (0), 
has been replaced in the new basis by a constant a(r ) _1 , no worse than the flat space or 
conformal wave function. Working in this modified vacuum, the infrared singularities will 
be removed, and quantities such as the propagator in our new vacuum will be well-defined. 
For example, the anticommutator G\ defined with respect to the state |r ) becomes 

G 1 (r,r',f-f')= /^kKM^+H" ( 28 ) 
For A; | To | > 1, the integrand is unchanged from the de Sitter- invariant vacuum, but for 



fc|T | < 1, the integrand is modified in such as way that it tends to a constant. If one considers 
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large spacelike separations, e.g., r = t', \x — x'\ — > oo, the behavior will be dominated by 
wave numbers k < l/\x — x'\. For this quantum amplitude, one is not restricted to the 
causal region, so one may consider the regime in which \x — x'\ ^> |to|. In this regime, the 
dominant contribution will come from the nonadiabatic modes where the wave functions 
may be approximated by Eq. (0). Thus, for large separations, G\ will fall to zero, 

G 1 {r,r,x-x')r, { R2 % 2 . (29) 
(x — x'Y 

Thus, this vacuum state satisfies cluster decomposition of the fields. On the other hand, 
the correlator G\ will increase logarithmically out to separations on the order of | To | before 
turning over, so that pairing only significantly changes the behavior outside the causally con- 
nected domain. Nevertheless, this is important for understanding the behavior of radiative 
corrections. This damping of long wavelength modes could also have implications for the 
sort of vacuum fluctuations that emerge at the completion of an inflationary epic. 

At the other extreme, consider the behavior of G\ at short distances. Of course, the 
correlator blows up as in flat space, but the first contribution to the short distance expansion 
that depends on the structure of spacetime comes from (<p(r, x) 2 ). By translation invariance 
of de Sitter space in conformal coordinates, we may choose x = 0. Of course, an operator 
such as <f) 2 requires renormalization in any case. However, it has been argued pO, pTL E2L ffl 



that, in the usual, de Sitter vacuum, the renormalized vacuum expectation value grows as 
In \t\ owing to the infrared divergence in this state. In our new vacuum \tq), this is not 
IR divergent. Thus, any log |t| dependence must then come from ultraviolet (conformal) 
momenta, so that it becomes a question of renormalization. Since <f> 2 is not an observable 
operator in the present theory, it is something of a moot point, but we will return to it in 
subsequent sections. 

One might entertain other states that are equally good candidates as vacuua. For exam- 
ple, in order to make the transition from the nonadiabatic to adiabatic modes smooth, one 
would be tempted simply to adopt the form f3 k given in Eqs. (p0|) and (^) for all k, not just 
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the adiabatic modes, but then to multiply the adiabatic modes by a damping factor, such 
as exp(— 6fc|ro|) for some constant b of order 1. Correspondingly, the magnitude of at would 
be modified so that the constraint Eq. QiTf ) i s satisfied. This seems to us to be a perfectly 
acceptable definition of an adiabatic vacuum, although probably more can be said about the 
transition between nonadiabatic and adiabatic modes. These issues would seem to depend 
on the details of the history of the de Sitter phase. 

The usual objections to defining the no-particle state by reference to diagonalizing the 
Hamiltonian do not apply to our construction. Normally, one would be concerned that there 
would be an infinite density of particles produced at all times later than the reference time 
ro, but that is not the case here since it is only the very low frequency modes that are 
significantly modified. Moreover, we have argued that one may not observe the nonadiabatic 
modes as particles; at best, their effects could be felt indirectly, for example, through their 
contribution to the energy density. In fact, we will show that their contributions to the 
energy density are small, vanishing with r as Tq 2 . 

The preceding considerations are naturally relevant to the vacuum state of quantum 
cosmological gravity (QCG), the theory of simply the gravitational field in the presence of 
a non-zero cosmological constant but without any matter. As mentioned earlier, the small 
fluctuations about the de Sitter background involve gravitons that include modes which, 
in the linear approximation, satisfy the same equation of motion as massless, minimally- 
coupled scalars.i Thus, we expect the vacuum state of QCG to resemble |r ), in which the 
non-adiabatic gravitons of opposite helicity and momenta are paired. With some additional 
effort, such a description can be made gauge- invariant. The situation is analogous to a BCS 
superconductor; the gap parameter is not gauge- invariant, so pairing would appear to break 
gauge invariance. But it can be shown that the Ward identities are fulfilled, so it is only the 
global [/(l)-symmetry that is broken. Similarly here, covariant conservation of the stress 
tensor is intact, even though global de Sitter symmetry is broken. Unlike the superconductor, 
8 This has been reviewed in Refs. @, ||, 0. 
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however, there are no Goldstone bosons associated with this symmetry breaking. In this 
respect, our situation resembles explicit breaking: even though the equations of motion 
remain de Sitter invariant, the action over the interval r > tq is not de Sitter invariant. 



3 Stress tensor in the state |ro,m) 

As another application of this vacuum, consider the renormalized, one-loop stress tensor 
T^y. Even in the de Sitter-invariant state |0, in), it is not beset by infrared divergences, and 
its vacuum expectation value is finite, although its value depends on the renormalization 
prescription. To gain insight into the effects of modifying the vacuum in this way, we evaluate 
its behavior at one-loop order in the state |r , in). A direct calculation of the unrenormalized 
contribution gives for the nonzero elements 

< r «> - i/(^K>i 2 -tI"H' (30) 

We rewrite the various terms in the integrand as, for example, 

K(r)| 2 = (l + 2\f3 k \ 2 ) |0 fc (r)| 2 -2^e{a fc ^(r)}, (31) 

where we used \ak\ 2 = 1 + \f3k\ 2 in the first term. The value when (3% = is of course 
the de Sitter invariant contribution, whereas the terms involving f3k represent the change 
arising from our modification of the vacuum. The de Sitter invariant contribution is com- 
pletely independent of r . It gives rise to quartic and quadratic ultraviolet divergences which, 
in a mass- independent renormalization scheme such as minimal subtraction, is completely 
removed by counterterms.! The process of renormalization also introduces a well-known 



anomalous contribution [24| to the stress-tensor which will be small for H 2 /m 2 Pl small. The 
remaining terms involving f3k depend on the initial time To and cannot be removed by renor- 
malization. As our interest is in the difference between the stress tensor in the de Sitter 



'People often adopt a renormalization scheme that retains a finite value for this piece, but it is arbitrary. 
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vacuum and the stress tensor our new vacuum |ro, in), we will ignore the anomalous con- 
tribution here. Thus, the change in the contributions to the one-loop renormalized stress 
tensor is 



1 



00/ 



2tt 2 Jo 



Vkol 



kdk { [\(3 k< p k \ 2 - 9£e (aMl)} + k 2 [\M k \ 2 - 3te (a k p k <fi)] } , (32) 



with a similar expression for (T&). Although the expressions for j3 k are quite complicated, 
we may use simple scaling arguments to obtain the general form of (T^). First of all, note 
that |fc0fc| 2 — H 2 (l + {kr) 2 ) and \<pk\ 2 — H 2 (kr) 2 . Since a k and f3 k are functions only of kro, 
we may rescale k — > k/r to see that (T 0Q ) is of the form A/t^ + Bt 2 /tq, where A and B 
are numerical constants that depend on the detailed form of the function fi k . Therefore, the 
energy density and pressure are 



(T °) =P 

-W> =p 



H 4 



+A 



B 



T 



A - 2 



T 



(not summed on 



(33) 
(34) 
0, pro- 



B r*" 

3 r 2 + 3 4_ 

It is easy to check that these forms obey the covariant conservation constraint T u 
viding an independent check on the relative signs appearing in the expression above. While 
our modified vacuum state contributes nonzero energy density (with (T^ u ) not proportional 
to the metric tensor g^ u ), these effects rapidly die away. (Recall that — Hr = exp(-Ht), so 
they decay exponentially in time.) The signs of the pressure provide insight into the two 
terms. The B-term is traceless and has the form of a relativistic gas that is diluted as a(r) -4 . 
The A-term has negative pressure and falls only as a(r) -2 . These are the amplified waves 
originally discussed by Grishchuk.0 Finally, as expected, as t — > — oo, these contributions 
vanish. 

Of course, because (T^) is non-zero, there will be feedback on the background metric, 
so the de Sitter form is no longer exact. However, for H 2 /m 2 Pl 1, this is a small distortion 
which, according to Eqs. (|33|) and (0), decays with time. 
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4 Modifications of the metric as an infrared cutoff. 



In this section we will pursue the same idea of introducing an infrared cutoff through the 
finite time of existence of de Sitter phase but in a different manner. Namely we will make 
an ad hoc assumption that metric has form (4) with the scale factor a(r) equal to 

a2 ^ = lk^ + lP (35) 
so that at large times, or r — > we recover de Sitter space. On the other hand assuming 
5 to be positive we see that as r — > — oo, or in the distant past, the metric becomes that 
of Minkowski spacetime. To be consistent with the Einstein's equations we tacitly assume 
a nontrivial matter distribution which would drive to the solution (^) (see also below). 
We will demonstrate that finite 5 does induce an infrared cut off both when one evaluates 
(O|0 2 |O) and (0|T M „|0). What might be most amusing, is that the infrared cutoff manifests 
itself in these two cases in different fashion. 

Since the metric has been modified we need to find the corresponding modes of the 
scalar field which would replace now solutions ([/]). We can try either to solve the equations 
of motion for the scalar field or to apply WKB expansion. The former approach has a 
shortcoming of being somewhat cumbersome while the latter, at first sight, is not suited at 
all to study infrared problems we are interested in. In fact the both ways work well. In 
this section we apply the high-frequency, or WKB expansion, while discussion of the exact 
solution is postponed until the next section. 

Let us first explain the general strategy of dealing with infrared problems within the 
WKB expansion and introduce to this end the necessary notation. We are looking for a 
solution 0^ with 3-momentum k which is represented as 

exp(zk-x) 

0k(x) = (2^a(r) Xk(T) (36) 



where Xk in massless case satisfies the equation (see, e.g., |J25 



^ + {k 2 - (£ - \)a 2 R} Xk = 0. (37) 
16 



Here R is the Ricci scalar which can be computed as 

R = a" 2 (3D + ^D 2 ^j , D = 2a/ a (38) 
and the dot denotes the derivative with respect to the conformal time r. The quasiclassical 



solution to (|37|) is 



Xk ^ (2W k (r)) 1 /2 ' (6J) 

Eq. ( f40|) allows to obtain H4 as an expansion in k~ 2 in a regular way. 
In particular, the first two terms in the expansion are well known: 

Wl = ^(l- ( ^f 2j W *)) (41) 

or 

k = (2tt)- 3 / 2 exp(zk • x - ikr) (r + 1 ^ - ^^ „.^ (42) 

while further terms can be obtained iteratively. 

What is special about de Sitter background and £ = 0, 1/6 is that these first two terms 
([42|) reproduce the exact result (see Eq. (|7|)). It is worth emphasizing that this termination 
of the expansion can be established within the quasiclassical expansion itself. The coefficient 
in front of k~ 3 term in the expansion (|42D does turn out to be zero if Ra 2 = 12/ V 2 . 

As a rule evaluation of a physical quantity involves integration over modes. Consider for 
example (O|0 2 |O). From Eq. (|42|) one immediately gets [25 



(0W0> 1 l^U-(A^E + o { ^)\. (43) 



(2tt) 3 7 2k \ 2k 2 

We will return later to the question of the ultraviolet regularization. In fact, the integrand is 
obtained in the high momentum limit, and the exact one should be infrared finite. However, 
the signal of the infrared cutoff may be inferred from the higher-order terms like k~ 4 . In case 
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that one starts with a quantity which is ultraviolet divergent (like (O|0 2 |O)), the following 
two simple rules seem to apply. First, all power-like ultraviolet divergencies can be removed 
through subtractions so that we have to concentrate on the log term. Clearly, the coefficient 
in front of the log is uniquely fixed by the k~ 2 expansion. Second, the next term in the 
expansion after the logarithmic one, i.e. the first term which diverges as a power in the 
infrared provides with an estimate of the infrared cutoff in the log. Indeed the infrared 
cutoff can be estimated by finding at which k 2 the term producing the log and the next one 
are of the same order. 

From this point of view, infrared problems in the de Sitter background arise not because 
subsequent terms in the k~ 2 expansion (see Eq. (f42|)) produce power- like infrared divergences 
but because for £ = 1/6 the high-frequency expansion terminates, and there are in fact 
no further terms in (f43|). As a result there is no infrared cutoff for the log and (O|0 2 |O) is 
undetermined. This is just the conclusion which was drawn from analyzing the exact solution 
for 4> k . Now we expect that by introducing 5^0 (see Eq. fl35[) ) we destroy the conspiracy 
of the coefficients which terminates the expansion ( pE2| ) and there emerges an infrared cutoff. 
We proceed now to check these expectations explicitly. To first order in 8 



Let us note that as far as Eq. fllPf ) is concerned such a modification of the metric is equivalent 
perturbatively to introducing positive £ ^ 0, 



Namely, negative £i? is like a negative mass-squared inducing an instability that leads to 
spontaneous breakdown of the de Sitter-invariant vacuum state. As a result, one cannot 
match the WKB expansion with the correct infrared behavior. If £ is positive, there is no 
danger of mismatch between the infrared and ultraviolet expansions. 




(44) 




35r 2 . 



(45) 



As is well known, flTTH there is an important difference between positive and negative £. 
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Keeping terms linear in 8 and expanding in k 1 leads to the following results: 

~H ~ (2tt)3 ^ ~ jfc + VP + PP + - 2jfc ~ fc2 + ifc3 + + >) ■ 

(46) 

Here we indicated the vanishing coefficients as well since they illustrate the statement above 
that the WKB expansion is perfectly consistent with the exact result for <pk (at 5 = 0). 

Now we are in position to apply the rules of defining the infrared cutoff mentioned above. 
Indeed, we see that the series in k~ 2 which represents (O|0 2 |O) does have further terms. The 
infrared cutoff, km can be estimated by equating two subsequent terms in the expansion: 

klk ~ 8k]%, k 2 IR ~ 5. (47) 



Note that the infrared cutoff (47) is time independent. Thus, for (O|0 2 |O) we find (omitting 
the subtraction term associated with the UV quadratic divergence) 

(O|0 2 |O) = l^lnA 2 ^, £=1/6. (48) 

According to the standard arguments the ultraviolet cutoff, Kuv does depend on time in the 
conformal coordinates, 

A uv ~ Ma(r) (49) 

where M is some mass parameter. Indeed, all the masses are scaled with a(r) in de Sitter 
background. If so, we reproduce the well-known result |20|, [21], [52] for (O|0 2 |O): 

(O|0 2 |O) ~ -^Inr. (50) 

We turn now to consideration of the vacuum expectation value (0 17)^0). The change 
brought by 5 ^ is a multiplicative factor (1 — 35 /2k 2 ) as far as energy and momentum 
associated with each mode 0k are concerned. If we turn to the integral over modes deter- 
mining the vacuum expectation value then this factor is crucial since the term proportional 
to 5 gives a logarithmically divergent integral over k. Keeping this log term we find 

(0|T 00 |0) = -3(0|r«|0> = -4^-ln^ (51) 

8ir 2 A IR 
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where Kuvi^-m are ultraviolet and infrared cutoffs and no summation over % is assumed. 
Moreover, as explained above, we do not keep power-like divergent terms. 

For the ultraviolet cutoff we still have fl4"5|). It is most interesting that the estimate ( pETD 
for the infrared cutoff is now changed. The reason is that the terms kept in the expansion 
([46|) are sufficient only to find the coefficient in front of the log term ([51]) while for the 
estimate of the infrared cutoff we need the next term in the k~ 2 expansion of the solutions 
0k- One can obtain the following estimate 

A 2 IR r 2 ~ 1. (52) 

It might worth noting that the time dependence exhibited by fl52p is so to say normal, 
since the WKB expansion is an expansion in inverse powers of kr. The fact that this time 
dependence did not manifest itself in (|4"?D is quite unique and is due to the termination of 
the WKB expansion in pure de Sitter background. 

Thus, we have for the contribution to the vacuum expectation value of the stress tensor: 

p = (0|T» cc -H*t% V = ~p. (53) 

Two remarks concerning this result are now in order. First, the vacuum expectation value 
is to satisfy the covariant conservation condition 

D M (0|2?|0) = 0, or r^-p = 3(p + p) (54) 

where the latter equation makes use of de Sitter metric, and the energy density p and pressure 
p are defined in Eqs. (|33[ , [34]). The use of de Sitter metric is justified since (0| X^,^ 1 0) starts 
from terms proportional to 5 and in the linear approximation one can neglect that 5^0 
otherwise. 

The condition ([54]) establishes a relation between time dependence and the coefficient of 
proportionality between p and p, in our case p = — l/3p. Moreover it is worth emphasizing 
that constraint ([54]) should be satisfied both for each mode with fixed k and for the sum over 

20 



the modes, or the vacuum expectation value. Since for a logarithmically divergent integral the 
relation between p and p for the sum over modes is the same as for each mode one immediately 
concludes that no extra time dependence can be associated with the \d.(Auv/^-ir)- This is 
a nontrivial check of general covariance which is passed by the estimates above. 

Second, Eq. ( p3[ ) implies that (0|T M1/ |0) is not proportional to g^ u which in de Sitter 
background would result in (0|T o|0) = — (0|Tjj|0). If we were considering pure de Sitter 
solution (i.e. 5 = 0) then g^ v would be the only tensor structure available, and Eq. ( |53D 
would not be possible. However, in the approximation linear in 6, we do have a new tensor 
structure emerging in R pv : 



By virtue of equations for the gravitational field the term proportional to 5 reflects the 
structure of T^ v of matter needed to produce solution (P5|). This new tensor structure is also 
manifested in Eq. (53). 

Thus, by assuming the scale factor to be given by (|35|) we have tacitly introduced non- 
trivial stress energy of matter which drives metric to such a solution. Then contribution 
(|53l) can be viewed as a radiative correction to the stress tensor of the matter. The plausible 
interpretation of this correction is that during the transition from Minkowski to de Sitter 
phases scalar particles are produced. 

In conclusion of this section, let us compare this method with that of the preceding 
sections involving the introduction of the initial time r . Note that, although technically the 
prescriptions for the infrared cutoff look very similar since they modify the situation in the 
distant past, the physics is actually quite different in these two cases. Indeed, in the latter 
case we assume that there is some matter distribution which causes (017)^0) ^ const ■ g^ u . 
In the former case, the deviations from pure de Sitter solution are entirely due to the initial 
state. 
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5 Exact solution for a modified metric. 



Here we will illustrate some of the results discussed above for the particular case of a modified 
de Sitter metric which was originally flat (at negative time infinity) and which admits an 
exact solution for the minimally coupled scalar field <fi. We chose the metric to be of the 
usual Robertson- Walker form with the scale factor: 

a (t) = 1 + exp(Ht) = [1 - exp(Hr)]- 1 (56) 

where the conformal time r is expressed through t as 

r = -H- 1 ln[l + exp(-Ht)] (57) 



As above the future infinity, t — > +oo, corresponds to r — > — 0. The scale factor ( |56"D 
cannot be created by any realistic source. For large positive time the source is close to the 
cosmological term and this can be achieved but in the past when both the energy and pressure 
density were close to zero the pressure should be negative and much larger by magnitude 
than energy. The latter cannot be realized by normal physical matter. Still this metric is 
convenient for the analysis of the quantum fluctuations in the de Sitter spacetime because 
the equation of motion (f|) can be exactly solved. The solution which behaves as exp(ifcr) 
when r — > — oo, has the form: 

<P k = y lK F(a,P,r,y) (58) 



where n = k/H, F(a, P,j;y) is the hypergeometric function as determined e.g. in ref. [p7 
a = — 1 + y/l — + ik, (5 = — 1 — \/l — k? + m, 7 = 1 + 2m, and y = exp(Hr). When r 
changes from — oo to —0, y changes from to 1. 

In the infrared limit when k — > and y is fixed, the solution is well defined, 

F -> 1 +m(-2y + y 2 ) (59) 

in contrast to expression ([?[) which diverges as 1/k. One would expect however that this 
divergence is recovered after a long de Sitter stage when y — > 1 (or r — > —0). This particular 
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hypergeometric function may be approximated near y = 1 as 



r(7) 



r(7-«)r( 7 -/3) 



2 + 7(y-i) + 2 ^ L ^(y-i 



\2 



(60) 



Once again the limit of small k is nonsingular and agrees with expression ( 59 ) for y — 1. We 
will observe the trace of the above mentioned infrared singularity in the ultraviolet expansion, 
k — > oo, of the solution. The result depends upon the way the limits k — > oo and y — *■ 1 are 
taken. For the case of y — *■ 1 first and then k — > oo, expression ( |60"D is convenient. We find 
from it (up to y-independent phase factor) that 



7T 3 + 2^ 
1 - -r- + 



(61) 



4k 24k 2 

The previously found solution for the exact infinitely long de Sitter stage is reproduced by 
the first two terms in the expansion in (y — 1) in the limit of large k. The infrared singularity 
at k = in the exact de Sitter solution (^) is evidently absent in the distorted metric and 
appears here only in the ultraviolet expansion. Note that due to the factor in front of the 
r.h.s. of Eq. flBTf ) the expansion of \F 2 \ contains not only even but also odd powers of 1/k. 
It would generate new types of divergences in perturbation theory. 

All these phenomena may be related to a singular character of the limiting transition 
from the metric ( |56|) to the exact de Sitter one. Indeed let us introduce a small parameter 
e (instead of 1 in Eq. (|5B])) which characterizes a deviation from the de Sitter metric at a 
fixed finite time, so the metric has the form: 

a(t) = e + exp(Ht) (62) 

By shifting the zero of time t = t' — ti so that exp(— Ht{) = e, we may rewrite a(t) as 
a(t') = e[l + exp(Ht')]. Now by rescaling the space coordinates dx = edx' we return to the 
old metric (|56"D. These transformations correspond to rescaling of the conformal momentum 
k — > k/e. For finite e this change of normalization is unobservable. But the de Sitter limit, 
when e — > 0, is singular and gives rise to the peculiar behavior mentioned above. 
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It may be also interesting to get ultraviolet expansion of the solution ( |58|) for a fixed y. 
To get this expansion it is convenient to present F in the form 

F = l-y+-f + ±g (63) 

K K 

The functions / and g can be found as expansions in inverse powers of k 2 : f = j\ + f 2 /n 2 + 
fs/n 4 ... (and the similar expression for g) by the perturbative (in 1/k 2 ) solution of the 
equations: 

(1 "</)/' + / = -i±» + ^ + ^/ (64) 

d-y) 9 ' + 9 = -i^/'-^/" (65) 

The terms which are logarithmically divergent in momentum are determined by the lowest 
order corrections f\ and g±. The latter are singular, though not too strongly, near y — 1 but 
this divergence is cancelled out, as we see in what follows, in the expression for (<f) 2 ). The 
calculation is straightforward and we get: 

h^-v-^-H^-v) (66) 

gi = V -- y -Hl-y)- l -{l-y)\n 2 {l-y) (67) 

The singularity in (1 — y) is getting stronger with the rising order in 1/k 2 . For example the 
most singular parts for f 2 and g 2 are respectively f 2 ng = 3/8(1— y) and g 2 ng = — 1/4(1 — y) 2 . 

Using these expressions we may calculate the vacuum expectation value of <fi 2 as a function 
of time. For a pure de Sitter universe, this quantity is known to rise as In r = t |2D|, ET1 



This result is ascribed to the infrared instability of a minimally coupled massless scalar 
field. We will calculate the difference (6(f) 2 ) = (0 2 (r)) — (<p 2 {— 00)). One may hope that this 
automatically takes care of infinities in flat spacetime. In the limit of small y one may write 
keeping only terms of order y 2 : 



(<t>\y)) 



d 3 k 



^2 2 V {2K 2 -l)y 



4/t 2 + l (An 2 + 1)(k 2 + 1) 
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(68) 



where y = exp(ifr). One has to subtract from this quantity the value at y = 0. The first 
term is quadratically divergent and to specify its value we have to define the regularization 
procedure. Usually it is just an arbitrary renormalization constant but in this case it becomes 
a function of time (or y). To avoid this time dependence one may change the variable of 
integration, k = a(r)p which means transition from the conformal momentum to the physical 
one. Written in this form the quadratically divergent integral in the expression for (f) 2 {y) is 
formally time independent and is cancelled out by the initial value in the flat spacetime. 

There are also logarithmically divergent terms which are absent in the flat case. For 
small y they can be found from Eq. (|68D which is exact with respect to the dependence on 
k and correspondingly is not singular at k = 0. Calculating the integral we find 

") = ^y( 1 + y)^K 2 max (69) 



Here n max is (unspecified) ultraviolet cutoff. For arbitrary y we may use Eqs. ( |63| , |66| , |67D 
giving expansion in inverse powers of k. Though the result is formally singular at k = 
we know that the effective infrared cutoff is K ~ 1. The logarithmic dependence on (1 — y) 



which we see in f 1 and g\ is cancelled out and we get exactly the same result (|B"9"|). The 



logarithmically divergent terms are explicitly y-dependent and this dependence cannot be 
cancelled by a redefinition of the ultraviolet cutoff as it was done above for the quadratically 
divergent part. If we assume that there is a physical cutoff realized, say, by some massive 
particles, we would expect the cutoff in conformal momentum to be proportional to the 
scale factor k = mo(r). This is because mass enters equations of motion in conformal 
coordinates always in the combination ma(t). In this case we get the result {5(f) 2 ) ~ y(l + 
y) ln(l — y). At the de Sitter stage when y — > 1, we get the well-known result (5(f) 2 ) ~ t 
plus a logarithmically divergent term. Note that the coefficient in front of the log is exactly 
R = Qy(l + y)H 2 . Expressed in terms of R the combination (1 — y) looks rather unusual: 



[1-V) = y/1 + 2R/3H 2 - 3. 

Let us turn now to the calculation of the energy-momentum tensor of the field in 
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background (1). Though it is a more UV divergent quantity than 2 , it has the nice property 
of covariant conservation which helps greatly in a proper definition of the divergent parts. 
Let us start from a more simple example when the background metric is flat in the infinite 
past and future and changes as a function of time in between. Assume that the scale factor 
a(— oo) = 1 and a(+oo) = aj. The mode decomposition of <fi has the standard form (8) with 
4>k(T~) = exp(-iuT) at negative time infinity and 4>kij) = [otk exp(— iur) + fa exp(iuT)f at 
positive time infinity. The Bogoliubov coefficients and fa satisfy the relation Eq. (|T8D . 
We will use the Heisenberg representation and will calculate the variation with time of the 
vacuum energy density. One has in conformal coordinates 



The physical energy density p is expressed through it as p = T u = T TT /a 2 . We will consider 
the expectation value (p(t)) averaged over the initial vacuum state. It is an infinite quartically 
divergent quantity, but the difference 5p = (p(+oo) — p(— oo)) is finite and is equal to 



It is the energy density of the particles produced by the time varying gravitational field. In 
this description we do not need a precise meaning of "particle" , the only essential quantity is 
the energy density. One can easily see that the pressure density of the final state is Sp = 5p/3. 
We use this simple example to illustrate that the finite part of the energy-momentum tensor 
averaged over the in-vacuum state is not proportional to the metric or to any other tensor 
associated with the problem and for the finite parts nonlocal effects are essential. Note also 
that the finite time difference 5p(t) = (p(t) — p(— oo) is quadratically divergent, and one 
needs curvature-dependent counterterms to eliminate these infinities. 

The situation is more complicated for the metric (|56"D when the spacetime is not flat in 
the infinite future. Still we may use the similar approach to calculate the evolution of the 
energy density of the field starting with the initial vacuum value which we assume to be 



Trr = \[{d T 4>) 2 + (d j( P) 2 ] 



(70) 




(71) 
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zero. The average value of the energy density over the initial vacuum state is 

1 r d 3 k u 



P 



a 2 (t)J (2tt) 3 2 



rv___ / T71 /\ , H y | (72) 



\ F 2 \ + ±^Sm(F*F') + -^-\F'\ 2 



k y ' k 2 

Here p = T tt is the energy density in the physical frame. Substituting expressions ( |63| , |66i|67D 
into this equation we get 

d 3 k uj 

— < i — a)- - - 

k 2 2k A 



P 



H 2 y 2 H\_ 

i , ^ 



:i - yf + -if - + ^rMy) 



a 2 J (2vr) 3 2 

The first two terms, quartically and quadratically divergent respectively, may be eliminated 
by renormalization (see discussion in Sec. 3). The renormalized quantities are at worst 
constants so their contribution into p is not rising with time. The last term is logarithmically 
divergent, but nonsingular at y = 1. The function K(y) might have in principle terms 
behaving like 1/(1 — y) or ln 2 (l — y), etc., as one can see from Eqs. ( |66| , |67l ) but all they 
have cancelled out. Hence the energy density of a massless minimally coupled scalar field in 
de Sitter background does not rise with time but goes down at least as I /a 2 and its back 
reaction on metric is not essential. Higher loop contributions probably will not change this 
conclusion, but this remains an open question. 

6 Conclusions. 

In this paper, we have provided a framework for resolving the infrared problems associated 
with the vacuum of de Sitter space. As previous workers have remarked, there is no solution 
without breaking de Sitter invariance. We have shown that the cosmologically well-motivated 
assumption that the de Sitter phase did not exist indefinitely is sufficient to remove these 
infrared divergences. The detailed modifications will of course depend on the particular 
history that precedes the era in which the vacuum energy dominates, but we have shown 
that the infrared problem can be resolved without reference to these details. 

We have offered several approaches that may be useful in different contexts. In the first, 
we suggested a treatment of the nonadiabatic modes that differs from the conventional one. 
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The single new parameter is the initial time tq before which the vacuum energy is presumed 
not to be dominant. In a manner analogous to the classic Lee-Nauenberg discussion of mass 
singularities, [|KS|] we argued that one must diagonalize the Hamiltonian before embarking 
on a perturbation expansion. We found this to be an economical formalism, in which the 
correlation functions differ significantly only at distances large compared to the distance 
that a signal might have traveled since the initial time tq. We showed that the correlation 
functions satisfy cluster decomposition, so that this is an acceptable framework for discussing 
quantum corrections. 

In the second approach, we explicitly modified the metric at early times, in such a way 
that, in the distant past, the space is Minkowski spacetime, and we explored the consequences 
at later times. While this cannot be solved analytically, a WKB solution provided insight 
into the nature of the infrared cutoff. In the third approach, we simply modified the scale 
factor in a simple manner, again so the past is Minkowskian, but in a way that the equations 
of motion can be solved analytically. 

Each approach has its strengths and weaknesses, but the first involving pairing has the 
advantage that it presumes little about the universe prior to the de Sitter phase, and the 
dynamical equations remain unchanged, so that their symmetry properties remain exact. 
Whether the transition between nonadiabatic and adiabatic modes can be refined in some 
general way, for example, through a path integral approach,0 is not clear. The pairing 
mechanism, as in superconductors, when applied to gravitons may have far-reaching con- 
sequences that we can not now foresee. Whether, in the case of gravity, this construction 
requires modification to take into account the self-interactions remains to be investigated. 
The approach in which the metric was explicitly modified may be useful when we know the 
universe is approximately Minkowski spacetime in the distant past. It also illustrates how 
the WKB method may provide insight into the nature of the IR cutoff in situations where 
the nonadiabatic modes cannot be so simply treated. 
10 See, e.g., Ref. §§. 
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We emphasize that the considerations of this paper do not resolve the issue of whether 
there might be long-time singularities that destabilize the de Sitter metric. || [7|, |8| It remains 
to be seen what consequences, if any, our results have on interacting field theories such as 
those considered by Ford|4[] in which the behavior of 2 directly influences the stress tensor. 
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